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Abstract 

We discuss photon emission from particles decelerlated by the cos- 
mic expansion. This can be interpretated as a kind of bremsstrahlung 
induced by the Universe geometry. In the high momentum limit its 
transition probability does not depend on detailed behavior of the 
expansion. This may play an important role when massless particle 
emission is discussed in the early Universe. 



1 Introduction 



Particles in the expanding Robertson- Walker Universe are decelerated in 
the comoving frame. Ignoring backreaction, they obey the geodesic equation 
and lose their physical momentum Pphys as 

PconJ 

Pphys = — ^ 0, (1) 

where Pconf is a conserved conformal momentum and a is a scale factor 
growing in time. In general particles in deceleration can emanate radiation, 
or some massless particles. We call this process geometric bremsstrahlungf] 
due to the cosmic expansion. 

Many analyses on phenomena in the early Universe have been performed 
so far using results of high energy particle physics and proposed a lot of inter- 
esting features of the Universe[0]. However they are based on calculations of 
transition matrices in the fiat spacetime, emphasizing the fact that rates of 
interactions are much larger than expansion rate given by Hubble parameter, 
and no careful attention seems to be paid to the geometric bremsstrahlung 
process, which actually gives no contribution in the fiat spacetime. 

We shall argue in this paper that geometric bremsstrahlung may play 
impotant roles in the early Universe. In the section 2, emission of electro- 
magnetic wave from a classical charged particle in the expanding Universe 
is discussed, taking account of backreaction. Our argument in the classi- 
cal level suggests that the emission rate may not be simply ignored and the 
damping time may nearly equal to expansion time. Thus in the section 3, we 
treat photon emission from charged a particle quantum mechanically. High 
momentum limit of the transition probability can be obtained analytically. 
We stress that massless limit should be treated carefully and is nontrivial. 

In this paper, we adopt the natural units, the light velocity c = 1 and the 
Planck constant ^ = 1. Signature of metric is taken as (+, — , — , — ). 



^ This process is also called by DeWitt and Brchme electro-gravitic bremsstrahlung in 
ref §. 
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2 Classical Geometric Bremsstrahlung in the 
Early Universe 



The radiation reaction has been neglected in the study of the early Uni- 
verse. However particles are deaccelerated due to the cosmic expansion and 
thus they will emit the radiation. If the damping time due to the radiation 
reaction is comparable to the expansion time, the effect of the radiation reac- 
tion may not be simply neglected. We shall study in detail this phenomenon 
in the case of classical charged particles. 

The study of the radiation reaction for a charged particle has a long 
history. The first relativistic calculation was performed by Diracp|. His cal- 
culation has been generalized by DeWitt-Brehme|]^ for the motion in gravita- 
tional field. They have shown that bremsstrahung induced by the spacetime 
curvature which we call geometric bremsstrahlung occured in addition to the 
usual radiation damping. The effect is nonlocal in general which is caused by 
the so-called tail term in the Green function. It was Hobbs[^] who corrected 
the result of De Witt-Brehme and pointed out that the tail term vanishes 
identically in the case of the conformally fiat spacetimes. His equation of 
motion for a particle with 4 velocity u^, mass m, charge e without external 
electromagnetic field may be written in the following form in conformally flat 
spacetime. 

where D / Dt is the absolute derivative along the worldline of the particle with 
T the proper time and the exp(2r2) is the conformal factor, g^i, = e^^rj^j^ with 
rj^y the fiat Minkowskii metric. 

Here we are interested in the radiation reaction induced by the cosmic 
expansion in the early Universe and thus we will restrict ourselves to the case 
where the conformal factor depends only on the time variable. We shall take 
a different approach from Hobbs and evaluate explicitly the damping time 
scale due to the geometric bremsstrahlung in the early Universe. 

We shall take the standard form for the action of a charged particle with 
mass m and charge e in the gravitational field. 
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S^-mj ^g^,x>^x-dT - e j A^x'^dr -^j d^x^g''^ g'^'' F^^Fp,. 

The dot denotes the derivative with respect to the proper time r. The 
equation of motion derived from the above action may be written as follows. 

Since we are interested in the early Universe, we may neglect the spatial 
curvature and thus take the spatially flat Robertson- Walker model as our 
background geometry, 

ds'^ = a^{ri){drf — dx"^) = dt^ — a^dx ^. 

Since the metric is conformally fiat, it is convenient to work in the con- 
formally related flat spacetime. Defining the conformally related proper time 
dTf — a'^dr, we shall define the conformally related 4 velocity 



dxi^ 



dTf 

Then the equation of motion may be written as follows. 

m = ma'^ \ -'-u^' 5'^ = ea -^rj^^FaBu'. 

Dt \dTf a dTf a J 

where the prime denotes the derivative with respect to the conformal time 
rj. By using the self field of the particle in the right hand side of the above 
equation, we shall obtain the radiation reaction force. 

Before calculating the reaction force explicitely, let us compare the timescale 
due to the radiation damping with that due to the cosmic expansion to see 
the importance of the radiation reaction in the early Universe. The damping 
time may be roughly evaluated as follows. 



r 



1 dEconf 



Econf adr] 




2 

26 I Pconf dd^ PphysH 



2 



Pconf 3a V / ^PconfO' \ "^^^ dr) J 3 m 
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where 

1 da 
H = — — 
a at 

is the Hubble parameter, Pconf is the conformal momentum and we have used 
the fact that the physical momentum Pphys = 'mu = Pconf/d decays as a~^ if 
the radiation reaction is neglected. Thus the ratio between the Hubble time 
texp = and the damping time is 



Pphys H . 



This ratio is much larger than the unity for a relativistic particle at sufficiently 
early times in the Universe. Thus the radiation reaction may not be simply 
ignored and might play an important role in the early Universe. 

For the calculation of the reaction force, we shall need the field equation 
derived from the above action, 



Taking the following non-covariant gauge 

we arrive at the field equation which has the same form with that in the fiat 
spacetime. 

Then the calculation by Dirac applies here and we obtain the standard 
expression for the reaction force in the fiat spacetime. 



d'^u^ ( du^"^ 



dTf \ dr- 



It can be shown by a direct calculation that our expression of the equation 
of motion with radiation reaction in the conformally related fiat spacetime 
coincides with eqn(0) when transformed back to the original physical frame. 
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In order to see the effect of the radiation reaction exphcitly, we shall 
forcus our attention to an 1-dimensional motion. Then the above equation 
is simplified as 



- — [amu) = -e 

CLTf 3 



d u u ( du 



dr/ 1 + ^2 \dTf 



Without the radiation reaction, the conformal momentum Pconf = clitlu is 
conserved as expected. 

Now we shall rewrite the above equation using the conformal momentum 
Pconf and the background time dt = adrj, 

d'^ Pconf / , 3m \ dpconf dH 

,^2 = \ H+ . = — — ^ + -7-Pconf- (3) 

dt' V 2e^^l + ipp,ys/my) dt dt 

Notice that there will be no classical geomotric bremsstrahlung in the case 
of de Sitter expansion, namely H = const. We shall be interested in the 
relativistic case in the early universe, namely 

Pphys > m. 

Then the second term in the coefficient of dpconf/dt in eqn(||) is negligible. 
Thus when the particle is relativistic, its evolution is governed by the reaction 
force only and the Hubble time will be the only available time scale in this 
situation. In fact, the solution in this case may be written as follows. 

Pconfit) =Po(l- H{to) dt' exp ^* H{x)d:^ ^ exp l^j' H{t')dt' 
where we have taken the following initial conditions; 

Pconfit = to) = po, -^-^(t = to) = 

The second condition expresses the fact that the reaction force is absent 
at the initial time. The solution shows that the momentum decays in the 
Hubble time. Thus this process should not be simply neglected. However 
the above conclusion is obtained as a classical effect and it is not clear if the 
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geometric bremsstrahlung is still effective when the quantum effect is talcen 
into account. We shall discuss quantum geometrobremsstrahlung in the next 
section. 

3 Quantum Geometric Bremsstrahlung 

As argued in the section 2, geometric bremsstrahlung may work effi- 
ciently classical mechanically in the early Universe. Whether this notable 
process survives or not after taking quantum effects into account is rather 
nontrivial and this question will be addressed next. 

To define well-behaved quantum amplitudes in the expanding Universe, 
we consider spacetimes with Minkowskian in- and out- regions. The way of 
expansion is chosen arbitrary. The scale factor is described as 

a{v) = C{Xv) (4) 

where t] is the conformal time, is a constant exhibiting typical time 
scale of the expansion. The function C{x) in eqn(|[) is arbitrary except the 
following constraints, 

C{x ~ -oo) = 6, (5) 
C(x~oo) = 1, (6) 
C{x) > 0, (7) 

where b is some positive constant describing ratio of initial scale factor to 
final one. 

Consider first the photon emission process in massive scalar QED with 
conformal coupling to the background curvature. The action reads 

S = J d^x {{Vf, + teAf,)^*{V^ -ieAf")^ 

+ {lR-m^)<^*<l>-]F^^'F^^. 
4 / 

The photon emission process is prohibited in the flat spacetime by energy- 
momentum conservation. However in the expanding spacetimes the energy 
conservation law gets broken and the transition can take place. The transi- 
tion amplitude is given in the lowest order of pertubation such that 

Amp = -ie ( d'^x^/^i (^^V''*, - V^$}$i) A*^ (8) 
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where is initial (final) mode function of massive charged scalar field 

and y4* is final mode function of electromagnetic field. The scalar mode 
functions satisfy 

^ ^ " (9) 



+ - -i? $ = 0. 
6 J 

Redefining the field $ = $ ■ a(r]\ the wave equation becomes the Klein- 
Gordon equation with a time-dependent mass, 

+ m^aijif) !• = 0. (10) 
Here we introduce g'^.ir]) and g^^'irf) satisfying a Schrodinger-type equation, 



drf- 



m a(r]y 



(11) 



with the boundary conditions in the asymptotic in and out regions as 

exp (^—irj^Jpi"^ + m?h'^ 



9fM 



Y^(27r)32y^2^m262 

exp {^—ir^^p^ + rr 



(27r)32^^ 



(77 00), 



(77 ~ 00). 



They also satisfy the following normalization condition. 

1 



[9p9 p-9 p9p) 



(27r)3' 



(12) 



where the prime denotes the derivative with respect to rj. Then the mode 
functions can be expressed as 

■' a ■' a ^' 
a a 

The electromagnetic final mode function satisfies the Maxwell equation 
in curved spacetime, 

W{V,Al - V.A;) = 0. (13) 



Notice that in 4- dimensional conformally flat spacetimes this eqn(|T3|) can be 
reduced into the same form in the flat spacetime, 



d^{d,Al - d^AD = 0. 
Therefore we get easily the final mode function 

exp (i\k\r] — ik ■ xj 



(14) 



(27r)32|A;| 



(15) 



where e* is a helicity factor. 

Using the rescaled field, the amplitude, eqn(|D, is rewritten as 

Amp = -ie J d^xi (I'^a'^l'i - S^l-^l-i) A*^. 



(16) 



Because the photon emission lasts only during the epoch of expansion, the 
concept of the probability per unit time is ambiguious. So we shall use the 
transition probability itself. The transition probability W can be obtained 
from the amplitude such that 



W 



E ^ / d%d'k\Amp\' 



h=L,R 



(17) 



where the summation is performed over the photon helicity and V is the 
conformal volume of the space, which is cancelled by the factor {27tY6{0) 
coming from the conformal momentum conservation in l^mpp . After the 
helicity summation, the explicit form of W is obtained as follows. 



W 



(27r)6e2 



d^k 



(27r)32|fc| 



Q , — * 

d pfS{k + pf -pi 



{Pf + Pi 



yypf yp, ypf yp. 



By virtue of the Wronskian relation; 

d / out* in' out*' in 



iPf 



Pi^)9: 



2\ „out* in 



(18) 
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the form of W is more simplified such that 



4 P.' - 



{kpi 



(19) 



To grasp the behavior of W in the \f)i\ oo hmit, we first argue a case 
with scale factor 

a(r7) =e(77) + 6e(-r/). (20) 
Then the exact mode functions are derived as 



9^: 



ei-v)- 



exp y—if]yPi + rn^h^ 



A{pi) exp ( -i-qJpi'^ + ) + B{pi) exp ( ir]Jf)^^ + iv? 
+e(r^) ^ ^ ^ ^, (21) 



y(27r)32y^+m262 



e(r^) 



exp 



-0(-^) 



^(271)32^^)2^^ 

y4(pj) exp {^—irj^'p^ + im?}p^ + B{pf) exp {irj^p^ + 7ii?lP^ 



y(27r)32y^7T^ 



(22) 



where 



5(p1 



^ 2 _|_ ^2 

Ip"^ + mP'W' 

p 2 -)- 77j2 



Substituting eqn (piD and eqn(|2^) into eqn ([l9D and taking the high momen- 
tum limit, IpjI oo, it is shown that the terms proportinal to B do not 
contribute to W because of the damping behavior of B. Notice that taking 
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the high momentum hmit, the energy conservation law almost restores in the 
following sense. 

\p^\-\k\ + \Pi-k\. (23) 

This can be read easily from the t] integration part in eqn ([T9| ) . Contribution 
from the region where eqn(|2^) does not hold is severely suppressed by the 
energy conservation factor. Using the polar coordinate decomposition pi-k = 
\pi\k cos 9 with k = \k\ and taking \pi\ much larger than m, it is easily derived 
that no contribution comes from the k integral region lying between |pj|and 
oo. Hence we get 

W{\pi\ ~ oo) 

,2 



2(27r)2 4-15 \p,\-kJo 



■1 -1 



k - Jpi'^ + m? + J{\pi\ - kf + 2\pi\k{l - cos 6) + w? 



■1 -1 



X 



k - ^Vi^ + + V - W + 2\p,\k{\ - cose) + 

(24) 

We need infra-red cutoff \pi\5 in eqn(p4D due to the existance of massless 
photon. This infra-red divergence is well known one in flat spacetime quan- 
tum field theories with massless particles and it should be cancelled by an 
infra-red divergence of the self energy term f^. The cutoff 5 is physically 
determined by resolving power of soft photon observation. The 9 integration 
in eqn(p4D can be straightforwardly calculated. After performing this inte- 
gration and taking the high momentum limit \pi\ oo, the k integration is 
simplified and we finally obtain 

w'(ip-.i--)=i:^('-4+*-i)(T^'4-2)- (^«' 



Note that we have taken the helicity sum in eqn(P3|). Instead, it is also 
possible to evaluate W{h) independently with a fixed photon helicity. For 
left and right handed helicity, each probability is the same, a half of W in 
eqn(p5|). 

Furthermore we can also obtain the analytic forms of W{h) in the spinor 
QED for the case of eqn(En). Because both of the charged fermion and photon 
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have degree of helicity freedom, 4 helicity contributions must be considered 
separately. The probabihty in the high momentum hmit for 1/2 hehcity 
fermion decaying into fermion with hehcitity hfermion and photon with helic- 
ity hphoton is denoted by W {1/2; hfermion, hphoton) and is given for each case 
as follows. 

■^"^^4-2|. (26, 



1^(1/2; 1/2,1) = 


8vr2 


W^(l/2; 1/2,-1) = 


'\ fin 


1^(1/2; -1/2,1) = 


87r2 1^ 


iy(l/2; -1/2,-1) 


= 0. 





1-62 52 

(28) 
(29) 

No infra-red cutoff 5 appears in eqn(^8D and eqn (|29|) because spinfiip of the 
fermion enables observers to distinguish the bremsstrahlung from the self 
energy process. 

There exists a very useful aspect of W in the high momentum limit. It is 
supposed that the results eqn(|25| ) ~eqn (|29|) are exact not only for the special 



way of the expansion given by eqn(|2^) but also arbitrary way satisfying 
eqn(|D ~ eqn(|^. This implies that Vr(|pj| 00) possesses a remarkable 
universality with respect to the ways of the cosmic expansion. This property 
may be explained as the Lorentz contraction effect from the view point of 
the high energy particle. Imagine a particle running in the comoving frame. 
Suppose that the Universe begins to expand when the particle passes through 
a point A and the Universe ceases to expand when the particle reaches point 
B. The particle catches energy from the expansion only while running from A 
to B. Taking the high momentum limit, the length between A and B contracts 
to zero in the rest frame of the particle. Therefore the particle cannot see 
the details of the way how the Universe expands and thus the universality of 
W crops up. 

To see the universality more quantitatively, we shall discuss the scalar 
QED with an adiabatically slow evolution of the scale factor a{ri) satisfy- 
ing eqn(^ ~ eqn(|^. In the zeroth order adiabatic approximation(WKB 
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approximation) the mode functions satisfying eqn(^) is written as 

exp 

9f; ^ 9f; ~ 



if) ■ X — i Jq dri'dp"^ + m?a{ri') 



^(27r)32y^ + m2a(r/)^ 



(30) 



Substituting eqn(|30D into eqn ([T9|) and introducing the polar coordinate de- 
composition; Pi ■ k = \pi\k cos 9, we get 

W ~ — ^ / dkk / rf^sin^^ 
2(27r)2 J\p^\5 Jo 



X 



oo exp 
dr] 

-oo 



ikr] — i P dri Jp^ + m'^a{ri'y + i P dri'u{pi — k, r]') 



uj{pi — k, ri)Jpi'^ + m'^a{riy 



(31) 

where uj{pi — k,ri) = {\pi\ — ky + 2\pi\k{l — cos 6) + rn^o? and \pi\5 is the 
inra-red cutoff. Consider the high momentum hmit in eqn (|3lD . As mentioned 
before nonvanishing contribution to W comes from the integral region where 
the momentum holds the relation eqn(^). Thus it is enough to restrict the 
momentum region between \pi\ and \pi\5. Here it should be searched which 
integral region of 6 contributes, accompanied with the influence of a{rj), to 
the nonvanishing value of W in eqn(|3TD. Due to eqn(|23|), only emittion to 
nearly forward direction {9 ~ 0) is permitted and especially the integral 
region of 9 satisfying 

< <0{m/\pi\) 
gives the scale factor dependence to the W . Several expansions like 



p^ + m?a{ri'Y ~ \pi\ + 



i\2 



yield finally 



W 



X 



dk- 



o 



d9 9^ 



2{2TTy J\p,\5 \pi\-kJo 
dr] exp » ' I 



\ 2\pi 



mk9' 



2{\p,\-k) 2{\p,,\-k)^ 



12 



=2 „l 



y f^^^^ s 

dyz / dzz 



X 



(27r)2 Js 1 — y Jo 
dfj exp 







1 - 1/ 1-1/ 



(32) 



where we change the integral variables in the following way, 

k = \Pi\y, 
m 

\Pi\ 

m\~ 

Note that the function C( ^^f'^ r/) in eqn(p2D approaches in the high momen- 
tum limit to a step function, 



C 



'2X\pi\^ 



fi] ~ Qif])+bQ{-f]). 



(33) 



Therefore the value of W for arbitrary adiabatical cosmic expansion satisfy- 
ing eqn(|p ~ eqn([^ must equal to the specified value for eqn(^D|), and the 
universality is surely realized. If we dismiss the adiabatic approximation. 



the mode functions have reflection wave terms like in eqn(^T]) and eqn (|22|) . 
However amplitude of the reflection waves vanishes in the high momentum 
limit and the universality are thought to survive. 

Here we have a comment on the rate of the geometric bremsstrahlung. 
Comparing with the classical results in the section 2, it is noticed from 
eqn(^)~ eqn(pp|) that interaction rate of quantum geometric bremsstrahlung 
does not so large compared with classical one. This is due to the fact that 
quantum effect smears position of the classical point particle and dilutes the 
charge density. 



Now it is worth considering implication of the results; eqn(l25| ) ^eqn (l29|) 
in connection with the conformal symmetry. Since the massless limit m — 
forces the speed of the particle to reach the light velocity, the universality 
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with respect to the way of the cosmic expansion is maintained. In the lowest 
pertubation of the QED, the massless hmit is shown to be equivalent with 
the oo limit. Therefore again the same results, eqn(|2S|)~ eqn(P5|), 

come up in the m — limit and W really possesses the non-vanishing value. 
One might naively expect for the massless case that the amplitude in the 
conformally flat spacetime vanishes as in the Minkowskian spacetime, by 
virtue of the conformal symmetry guaranteed at least classical mechanically. 
However this is not true unless 6 = 1 as argued above. 

Next let us discuss the rate of the geometric bremsstrahlung. Taking large 
expansion limit 6 ~ 0, the transition probablity can be typically expressed 
as 



0) 



N*e'^ 1 
0(l)^ln^. 

^ ' Att^ b' 



(34) 



where N* is the number of final modes. In our spinor QED model we take 
A^* = 2 due to contributions of eqns (p6D and (pTj). However in the standard 
model we have more particles and A^* ~ 0(10). Moreover extended theories 
(like GUT and SUSY) can gives us A^* ~ O(IOO). 

The process is thought to occur when W ^ 1. Thus when the Universe 
expands enough to satisfy 



ai 



exp 



47r^O(l) 



(35) 



the event will take place. If we specify a model of the Unvierse evolution, 
we can get the rate itself. For example let us assume that the expansion is 
dominated by the radiation, a(t) oc t^^^. Then we get the following estimation 
for the transition rate. 



1 



ti 



exp 



^71^0(1) 



~ e 



O(IOO) 



Hi, 



(36) 



where Hi = 1/ti is the Hubble parameter at the initial time. Thus we cannot 
neglect the quantum bremsstrahlung naively when A^* ~ O(IOO). It might 
be also worth reminding that rates of ohther gauge interection processes 
different from the geometric bremsstrahlung is much smaller than the Hubble 
parameter when the temparature of the Universe is higher than 0(10^^) GeV. 
Thus the bremsstrahlung may be the most dominant process in such a early 
era. 
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